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A new mechanism for the pion production in high energy reactions is suggested. It is related to 
a possibility for the direct production of the pions induced by instantons, topologically nontrivial 
glnonic excitations of the QCD vacnum. This mechanism does not require any fragmentation func¬ 
tions for the production of pseudoscalar mesons in high energy reactions with hadrons. We calculate 
the contribution of the new mechanism to the inclusive 7r°-meson production in high energy proton- 
proton collisions. It is shown that it gives the dominant contribution to the inclusive cross-sections 
in the few GeV region for the transverse momentum of the final pion with large rapidity. We discuss 
the possible applications of the new mechanism to the phenomenon of large spin effects observed in 
numerous high energy reactions and to particle productions in the relativistic heavy ion collisions. 


I. INTRODUCTION 

The inclusive production of mesons in high energy re¬ 
actions is one of the powerful tools to investigate the 
structure of strong interaction. At the very large trans¬ 
verse momentum the leading behavior of the inclusive 
cross-section should be dominated by the perturbative 
t-channel one-gluon exchange which leads to the well- 
known dependency da ^ ^/Pt- However, one cannot be¬ 
lieve in the validity of the perturbative QCD (pQCD) ap¬ 
proach in the small transverse momentum region where 
the nonperturbative QCD effects can play a crucial role. 
It is clear that the value of the transfer momentum for 
the applicability of the pQCD depends on the nonper¬ 
turbative dynamics of QCD. This dynamics is deeply re¬ 
lated to the complicated structure of the QCD vacuum. 
One of the powerful models to calculate the nonpertur¬ 
bative QCD effects in hadron physics, which are induced 
by the nontrivial topological structure of the QCD vac¬ 
uum, is the instanton model (see reviews The in¬ 

stantons describe the subbarrier transitions between the 
classical QCD vacua with the different topological charge. 
The existence of instantons is very important for hadron 
physics. For example, they provide a natural mechanism 
of the spontaneous chiral symmetry breaking (SCSB) in 
the strong interaction. As the result, large dynamical 
quark masses arise. One of the places where the SCSB 
effects might be important is the high energy reactions 
with hadrons. In particular, it was demonstrated in Q 
that in the few GeV region for the momentum trans¬ 
fer the instanton effects give a significant contribution 
to the high energy quark-quark scattering cross section. 
These effects come from the anomalous chromomagnetic 


quark-gluon interaction induced by instantons Q. Fur¬ 
thermore, this interaction leads to the quark chirality flip 
and might give an important contribution to the spin- 
dependent cross sections [1-Q. The problem with the 
pQCD description of the high energy inclusive pion pro¬ 
duction was discussed for the first time in It was 
mentioned that the cross section is well described by the 
pQCD only in the region of small xp for the fixed tar¬ 
get experiments. In jlflj . the large partonic intrinsic mo¬ 
mentum in hadrons was introduced to describe the data. 
However, it is not easy to justify the existence of such a 
large intrinsic partonic momentum from the confinement 
dynamics. The possible violation of the pQCD factoriza¬ 
tion in inclusive production of hadrons induced by high 
twist contributions was discussed by Brodsky with col¬ 
laborators (see [mil and references therein), but the 
microscopic mechanism of such a violation was not pre¬ 
sented. 

In this Letter, we suggest a new mechanism of the inclu¬ 
sive production of pions in high energy reactions. We 
call it the anomalous pion production (APP) because 
the formation of the pions happens at short distances 
due to instantons which have a much smaller size in 
comparison with the confinement scale. Furthermore, in 
our approach it is not needed to include in the calcula¬ 
tion any fragmentation functions which are related to the 
hadronization and, therefore, to the confinement dynam¬ 
ics. As the result, the value of the inclusive pion cross 
section is determined by the structure of the short range 
fluctuations of gluonic fields in the QCD vacuum. This 
mechanism breaks the pQCD factorization and might be 
a corner stone of various phenomena observed in high 
energy reactions in the few GeV range for the transfer 
momentum. 
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II. ANOMALOUS PRODUCTION OF THE 
PIONS INDUCED BY CHROMOMAGNETIC 
VERTEX 


It was shown that instantons generate a new type of 
the quark-gluon chromomagnetic interaction @ 




( 1 ) 


where /ia is the anomalous quark chromomagnetic mo¬ 
ment (AQCM), Mq is the dynamical quark mass, gs 
is the strong coupling constant, and is the gluon 
field strength. Within the instanton model the value of 
AQCM is (see for details 0,0) 


/^a — 


Z-n{MqPc) 

4as(pc) 


( 2 ) 


Therefore, the value of AQCM is determined by the dy¬ 
namical mass of the quark in the instanton vacuum. For 
example, for the dynamical quark mass Mq = 170 MeV 
from the mean field approximation [l[ and as{pc) ~ 0.5 
at an average size of instantons in the QCD vacuum 
Pc = 1/600 MeV“^ 0, we obtain 


Pa = -0.378, 


(3) 


which is very large in comparison with the Schwinger- 
type of the pQCD contribution to the AQCM 

~ -1-3 X 10-2. (4) 

It is evident that Lagrangian, Eqll] violates chiral sym¬ 
metry. Therefore, in |2l.ll4| the generalization of Eq|T]was 
suggested by inclusion of the pion field into consideration 
to preserve the explicit chiral invariance. The modified 
Lagrangian is 

Cl = (5) 

where = 93 MeV. The expansion of this Lagrangian 




FIG. 1: The diagram a) presents the general quark-gluon ver¬ 
tex generated by instantons for Nf = 2 case , b) corresponds 
to the case with one of the quark lines connected through the 
quark condensate, and c) describes the direct production of 
the pion from the instanton. 


Within the instanton model the two terms in this equa¬ 
tion can be represented by the diagrams b), c) in Fig.l. 
The diagrams which give the contribution to the anoma¬ 
lous production of the pion on the parton level are shown 
in Fig.2. Using the Sudakov parametrization of the four- 


k ^ 7T 



a) b) 

FIG. 2: The pion production induced by instanton a) in 
quark-quark scattering and b)in quark-gluon scattering. 


momenta of particles [T^ 

Pi = a,P2 -f ZiPi + Pi t, p'^ t ■ Pi,2 = 0, pft = -^,2 


we obtain the cross section of the 7r° production in the 
quark-quark scattering by incoming u- or d-quark 




9^0 


dz(Pk 

Pt [ 

™277r J 


a2(Pc) (-f+m2)2 3 


{pc\q\)cPq, (7) 


where F„{i) = — 2K2{t) is the form factor of the 

instanton [l3|, 2 is a fraction of the initial quark mo¬ 
mentum carried by the final pion, rug is the dynamical 
gluon mass related to the infrared behavior of the gluon 
propagator, and 


PirOqq j-, ■ ( 8 ) 

For the quark-gluon scattering (the diagram in Fig.2b), 
we obtained the cross section which is larger than quark- 
quark scattering case by a factor of 9/4. The instanton 
corresponds to the subbarrier transition between vacua 
with the different topological charge. Therefore, the sin¬ 
gle instanton approximation, which we are using, is cor¬ 
rect when the invariant mass of the partonic system pro¬ 
duced by the instanton does not exceed the height of 
the potential barrier between these vacua. This height is 
given by the energy of the so-called sphaleron Esph = 
37r/(4as(p)p) (see, for example, 0)- For \/pc = 0.6 
GeV and a{pc) = 0.5 we obtain E^ph = 2.83 GeV, 
which is rather a large value. By using the condition 
{k+p'iY = Mjr < E'^pf^ and the relation 


up to the first order in the pion field gives 


28| 


{zq-kf=z{l-z)M‘^, (9) 


Cl 


■ 9 spa 
'^q 


qa^'^eqC^ 


9 s Pa 
iMqE^ 


( 6 ) 


where k is the transverse momentum of the pion, we fi¬ 
nally obtain the APP cross section in the quark-quark 
scattering 
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dz(Pk 


n2TT 


9Tr°qq 2 


(Pc)2« 


dy / d(p- 


where 


F{x, y, z,(p) = z{l — z)y + + 2a;-y/y-^/z^T^^^cos ip 

( 11 ) 

and X = pc\k\, y = This cross section can be 

compared with the leading order pQCD cross section 


zd(jm^'}'}^° 

dzd^k 


{pQCD) = 


Dq\z) da 
TTZ dq^ ’ 


( 12 ) 


where D'^ (z) is the fragmentation function and 
da/dq^ = 87ra^(g^)/(g^ — TOg)^/9. We should mention 
that in the kinematic region of few GeV for the trans¬ 
verse momentum of the pion and at high energy one can 
neglect the change of the longitudinal momentum of the 
quark in partonic subprocess in the pQCD. In this case, 


(x, y, p)F{x, y, z, ip)F^{y'F{x, y, z, ip)/z) 
{F{x,y,z,(p) +z'^mlplY 


z(l - z), 


( 10 ) 


the meaning of z becomes the same for both pQCD and 
APP cases. In Fig.3, the result of calculation of the cross 
section of the APP, EplTOl in the quark-quark scattering 
is presented [s^- We compare it with the pQCD cross 
section, EqUll calculated with the strong coupling con¬ 
stant given by the analytical pQCD 


= IT 


1 A2 

-k 


Po[log{QyA^) A2-Q2J 


(13) 


where = -q^, /3o = (33 - 21V/)/127r, A = 250 MeV, 
and Nf = 3 were used. The fragmentation function 
D^° (z) was taken from [13, set KKP-1. The dynami¬ 
cal gluon mass nig ~ 0.65 GeV was fixed according to 
the result of the lattice data in Eliall. 
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FIG. 3: The z dependency of the pQCD (black line) and the APP (blue line) cross sections, in the units phjGeV^^ for 
the different values of the pion transverse momentum: the left panel kt = 1 GeV, the central panel fct = 2 GeV, and the right 
panel kt = 3 GeV. Here, the relation kt = |fc| « \/—q^ was used. 


Let us discuss, in the beginning on a qualitative level, 
the features of kinematics of the pQCD and the APP 
contributions to tt^ production in the fixed target ex¬ 
periments for a few GeV transverse momentum of the 
final pion where the big difference from the pQCD pre¬ 
diction was observed. In this case, the main contribution 
comes from the scattering of the valence quarks which 
carry about < x >~ 0.2 momentum of the initial proton. 
Therefore, it follows from Fig.3 that at A:* = 1 GeV the 
APP contribution starts to be dominant at xp ~ 0.08 
[ 13 , at kt = 2 GeV the APP contribution dominates at 
xp ~ 0.11, and at kt = 3 GeV the APP above xp ~ 0.15 


gives the main contribution to the cross section. So we 
come to the conclusion that the pQCD might give the 
main contribution only in the region of very small xp. 
Furthermore, the APP contribution is very big even at 
large kt in the large xp region. To calculate the APP 
cross section at the hadron level, we use an approach 
similar to the pQCD approach (see [13 )■ In particular, 
in the high energy limit we will neglect all masses of the 
hadrons and partons in the kinematics of the inclusive 
pion production. In this way, for the APP contribution 
to the inclusive cross section of the production in the 
proton-proton scattering, we have 
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E^da 

d^k 



^A—^a(,^a7 d B—^b(,^b: d ) 


Zdcr“^ , 

- -{Z 

dzd'^k 


- ^C7 kt)^ 


(14) 


where GA,B^a,b are parton distribution functions 
(PDFs) and the sum over the different types of par- 
tons is carried out. Additionally, we should use the scale 
/j,^ Ri 1/p^ in our nonperturbative calculation. 




FIG. 4: The APP cross sections for inclusive 7r° in comparison 
with the ISR (CERN Intersecting Storage Rings) fixed target 
and the STAR data for kt ~ 1 T 3 GeV at large 
pseudorapidity. 


In Eg fin 

a;™” = a;i/(l-a;2), = X2Xa/{xa - Xi), 

Xi = -XT/tan{0cml‘^), X2 = -XTtan{ecm.l‘^), 

Zc = X2lXh + XijXa, (15) 

where Xt = ^ktj^fs and 9cm is the c.m. angle of the 
produced pion and the GRV98 NLO parametrization at 
^2 = 0.4 GeV^ was used [T^. This value of the scale is 
very close to the scale of our instanton based calculation 
1/p^c — GeV^. The results are represented in Fig.4 in 

comparison with the experimental data. As one can see, 
the APP mechanism provides a rather good description 
of both fixed target and collider data at xp > 0.2 and 
few GeV value for the transverse momentum of the pion. 
Some deviation from STAR data might be related to the 
PDF uncertainties. On the other hand, it is known that 
in this kinematic region for the fixed target experiments 
the pQGD contribution is very small and cannot explain 
the data without introducing large intrinsic parton mo¬ 
mentum I, 111 in the leading order calculation. The 
NLO pQGD predictions also fail to explain these data. 
We would like to emphasize that it is very difficult to 
uncerstand in the conventional approach why the pQCD 
works well at RHIC energy but fails to explain large xp 
ISR data. In both the cases the energy is much larger 
than any other scales, i.e., hadron masses, transverse mo¬ 
mentum of partons and hadrons, transfer momentum in a 
partonic subprocess, etc. Furthermore, the contribution 
from gluonic gg -A gg, which could be one of the candi¬ 
dates to explain this difference, is expected to be small in 
the forward rapidity region even for the RHIG energies 
[ 2 ^. It might be that a rather good agreement for the 
collider case [IJ is related to the internal uncertainties in 
the pQCD calculation due to uncertainties in the value 
of factorization scale and in the shape of fragmentation 
functions. Indeed, the comparison with the STAR data 
shows rather large sensitivity of the NLO pQGD predic¬ 
tion to the choice of the fragmentation function and 
to the factorization scale . Therefore, we believe that 
the complete picture of the inclusive particle production 
in the strong interaction should include both pQCD and 
APP mechanisms. 

III. CONCLUSION 

We discussed the new mechanism for the inclusive pro¬ 
duction of the pions in the hadron-hadron interaction. It 
is related to the anomalous quark-gluon-pion coupling 
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induced by the nontrivial topological structure of the 
QCD vacuum. It is shown that this APP mechanism 
gives the dominant contribution to the inclusive pion 
production for the xp ^ 0.1 in the few GeV region for 
the transverse momentum of the pion. Furthermore, the 
APP violates factorization in the inclusive pion produc¬ 
tion and allows one to calculate the cross section with¬ 
out any fragmentation functions. The APP mechanism 
could be extended to the other pseudoscalar mesons. In 
the case of K- and 77 -meson production, it is needed to 
introduce the suppression for the chromomagnetic inter¬ 
action due to the large mass of the s-quark. Addition¬ 
ally, one should include the effect from the 77 — 77 ' mixing 
in the consideration of 77 -meson production. It is evi¬ 
dent that the APP mechanism should give the contri¬ 
bution to other observables in the pseudoscalar meson 
production in high energy reactions with hadrons. We 
should emphasize that the anomalous quark-gluon-pion 
coupling, EqlSl flips the quark spin. Therefore, it should 
contribute to the large single spin asymmetries (SSA) in 
the inclusive pseudoscalar meson production observed in 
the different high energy reactions with hadrons. The 
other longstanding problem is the large transverse polar¬ 
ization of A polarization observed in different high energy 
reactions (see discussion in In our approach, 

such polarization might come from quark spin-flip in the 
A'-meson production by the APP-type mechanism, e.g. 
u + q{g) —>■ K'^ -I- s -I- g(g), and following fragmentation 
of the strange quark to A. In this case, the SSA in the 
A'-meson production and the transverse polarization of 
A should be deeply related to each other. However, the 
calculation of SSA is a much more complicated task due 
to the loop diagrams contribution even in the leading or¬ 
der. The APP mechanism is completely different from 
the mechanism related to the first term in EqI5] which 
was considered in ii at the partonic level. However, 
at the hadron level one should include the fragmentation 
function into consideration. It is clear that such a contri¬ 
bution to the pseudoscalar meson inclusive cross section 


should have approximately the same xp dependency as 
the pQCD contribution which cannot solve the problem 
of the pion production at large xp. Our separate task 
is to study the effects of the APP mechanism to the in¬ 
clusive production of the pseudoscalar mesons in high 
energy heavy ion collision, in particular, to the so-called 
nuclear modification factor Raa- In the case of pQCD 
production, the size of the pion is fixed by the fragmen¬ 
tation dynamics and should be rather large, ~ 1 

fm. For the APP mechanism the size of the pion is small, 
rApp ^ ^ Q 3 the final state in¬ 

teraction cross section of the pion in the nuclear matter 
for the APP process is about one order of magnitude 
smaller in comparison with the pQCD case. Therefore, 
the pQCD and the APP mechanisms should have dif¬ 
ferent nuclear dependency. This difference, for example, 
might be responsible for the peak in Raa observed at 
RHIC and LHC in the few GeV region for transverse 
momentum of the final mesons. The work in these direc¬ 
tions is now in progress. 


IV. ACKNOWLEDGMENTS 

We would like to thank A.E. Dorokhov and S.B. 
Gerasimov for useful discussions. We are very grateful 
to Umberto D’Alesio for providing us with the tables of 
experimental data for the cross section of the inclusive 
7 r° production in different high energy reactions. This 
work was partially supported by the National Natural 
Science Foundation of China (Grant No. 11035006 and 
11175215), and by the Chinese Academy of Sciences visit¬ 
ing professorship for senior international scientists (Grant 
No. 2013T2J0011). This research was also supported in 
part by the Basic Science Research Program through the 
National Research Foundation of Korea(NRF) funded by 
the Ministry of Education(2013RlAlA2009695)(HJL). 


[1] T. Schafer and E.V. Shuryak, Rev. Mod. Phys. 70 (1998) 
1323. 

[2] D. Diakonov, Prog. Part. Nucl. Phys. 51 (2003) 173. 

[3] N. Kochelev, JETP Lett. 83 (2006) 527. 

[4] N. I. Kochelev, Phys. Lett. B426 (1998) 149. 

[5] N. I. Kochelev, JETP Lett. 72 481 (2000) 

[6] N. Kochelev and N. Korchagin, Phys. Lett. B 729, 117 
(2014) 

[7] D. Ostrovsky and E. Shuryak, Phys. Rev. D 71 (2005) 
014037 hep-ph/0409253 . 

[8] Y. Qian and I. Zahed, Phys. Rev. D 86 (2012) 014033 

[9] C. Bourrely and J. Soffer, Eur. Phys. J. C 36, 371 (2004) 

[10] U. D’Alesio and F. Murgia, Phys. Rev. D 70, 074009 
(2004) 

[11] F. Arleo, S. J. Brodsky, D. S. Hwang and A. M. Sickles, 
Phys. Rev. Lett. 105, 062002 (2010) 

[12] S. Brodsky, G. de Teramond and M. Karliner, Ann. Rev. 


Nucl. Part. Sci. 62, 1 (2012) 

[13] N. Kochelev, Phys. Part. Nucl. Lett. 7, 326 (2010) 

14 J. Balia, M. V. Polyakov and C. Weiss, Nucl. Phys. B 
510, 327 (1998) 

[15] V. N. Baler, E. A. Kuraev, V. S. Fadin and V. A. Khoze, 
Phys. Rept. 78 (1981) 293. 

[16] D. V. Shirkov and I. L. Solovtsov, Phys. Rev. Lett. 79, 
1209 (1997) 

[17] B. A. Kniehl, G. Kramer and B. Potter, Nucl. Phys. B 
582, 514 (2000) 

[18] E. Ruiz Arriola, P. O. Bowman and W. Broniowski, Phys. 
Rev. D 70 (2004) 097505 

[19] M. Gluck, E. Reya and A. Vogt, Eur. Phys. J. C 5, 461 
(1998) 

[20] Richard D. Field, ’’Applications of Perturbative 
QGD”,Addison-Wesley Publishing Gompany, Frontiers 
in physics, v77 1989 



6 


[21] D. Lloyd Owen, G. W. Abshire, G. Finocchiaro, 
P. D. Grannis, H. Jostlein, R. D. Kephart, R. Thun and 
G. Bellettini et al, Phys. Rev. Lett. 45, 89 (1980). 

[22] J. Adams et al. [STAR Collaboration], Phys. Rev. Lett. 
92, 171801 (2004); Phys. Rev. Lett. 97, 152302 (2006) 

[23] S. Kretzer, Acta Phys. Polon. B 36, 179 (2005) 

]24] M. Boglione, U. D’Alesio and F. Murgia, Phys. Rev. D 
77, 051502 (2008) 

[25] I. Arsene et al. [BRAHMS Collaboration], Phys. Rev. 
Lett. 98, 252001 (2007) 

[26] M. Anselmino, D. Boer, U. D’Alesio and F. Murgia, Phys. 
Rev. D 63, 054029 (2001) 

[27] M. Anselmino, D. Boer, U. D’Alesio and F. Murgia, Int. 
J. Mod. Phys. A 18, 1237 (2003) 

[28] We will neglect high order terms. Their contribution to 
the cross section is expected to be suppressed in a large 
Nc limit by a factor of 1/Nc, because TLr ~ y/Nc. Ad¬ 
ditionally, due to increasing of the number of the hnal 
particles, it should be suppressed at large xp in compar¬ 


ison with a leading contribution. 

[29] The production cross section of (tt”) in forward u- 
(d-) quark scattering is twice larger. 

[30] The pion production by our nonperturbative mechanism 
in the quark-gluon scattering is scaled by the same color 
factor 9/4 as in the pQCD case. Therefore, the inclu¬ 
sion into consideration of the quark-gluon scattering does 
not influence our conclusion. We also don’t consider the 
pQCD contribution coming from the gg gg subpro¬ 
cess, because it gives the sizable contribution only at very 
small value of z. The pQCD contribution from the glu- 
onic subprocess gg —>■ gg is expected to be small in the 
large xf region even for the RHIC energies [^ . 

[31] We do not consider the possible effect coming from the 
running dynamical gluon mass. 

[32] The Xf is part of the initial proton momentum carried 
by the final pion in the center of mass system in the 
longitudinal direction. 


